
restore,'./irdc_dist_model/bgps_rb3.sav',/ver


dI = rb3.del_Imir
imin = rb3.i_min
sigI = 0.3 

rad = dI*dI + sigI*sigI/imin/imin

value = dI*dI / rad

plothist,value,bin=0.01,color=cgColor('Black'),/ylog

in90 = WHERE(value GE 0.9, n90)
print,float(n90) / float(n_elements(rb3))

END
